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Abstract 

NDE (Near-dissociation expansion) including LeRoy-Bernstein for- 
mulas are improved by taking into account the multipole expansion 
coefficients and the non asymptotic part of the potential curve. Apply- 
ing these new simple analytical formulas to photoassociation spectra 
of cold alkali atoms, we improve the determination of the asymptotic 
coefficient, reaching a 1% accuracy, for long-range relativistic potential 
curve of diatomic molecules. 

1 Introduction 

The interaction between two distant atoms has been first studied by 
Van der Waals and London (for review see [Q). This topic is often 
discussed as a limiting case between the Hund case (a) and (c) 
13]. The study of such excitation transfers [U E] are related to long- 
range molecular states, where the electronical potential V(R) is fully 
described by the asymptotic coefficients 1 

V{R)^D + Y% (1) 

k 

1 LeRoy and Bernstein and some other authors use V(R) ~ D — J2 k Qfc. 



for sufficiently large internuclear distance R 0. Among these long- 
range states Stwalley et al. discovered very particular molecu- 
lar states: the so-called "pure long-range state", where both classi- 
cal turning points are in this asymptotic area. Great efforts have 
been devoted to the precise calculations of the asymptotic coefficients 



Semiclassical formulas in diatomic molecular spectroscopy are pow- 
erful tools (for a brief review see jlUUllj ). Several molecular properties 
as rotational or vibrational progression and kinetic energy are strongly 
determinate by the leading terms in equation ([T]): 



where m > n. In the course of this, article we shall suppose n > 2. 
In 1970, LeRoy and Bernstein pioneer work, based on the Bohr 
quantization formula, made possible to extract the leading coefficient 
C n from experimental vibrational progression. The LeRoy-Bernstein 
formula links the energy E of the vibrational quantum number v with 
the asymptotic behavior D + C n /R n of the potential curve : 



where \x is the reduced mass of the system and vd is the non- integer 
value of v at the dissociation energy D. This kind of near-dissociation 
expansion (NDE) semi-classical formula was extended to the rota- 
tional progression ^3] and the kinetic energy [T3j. The technique was 
also improved to include other coefficients in the asymptotic devel- 
opment ^H] and a quasi-complete NDE theory was established [TH] . 
The subject is still in progress: links with Quantum Defect Theory, 
scaling law for the density probability of presence of the vibrational 
wavefunction ^2 El > and extension to two coupled channels and Lu- 
Fano plots have been successfully investigated ^ElEIl- The S oa l °f 
this article is to improve part of the NDE theory. 

Experimental studies of the long-range states [20] have recently 
been renewed by the photoassociation (PA) spectroscopy of trapped 
cold atomic samples. Trapping and cooling of neutral atomic sam- 
ples, based on radiation pressure, are well established |U techniques 
that led to further spectroscopic developments. For instance, in 1987, 
Thorsheim et al. [221 proposed a new spectroscopic technique: the 
photoassociation process where a pair of free cold atoms absorbs reso- 
nantly one photon and produces an excited molecule in a well-defined 
ro-vibrational level. The first experiments were realized in 1993 in 
sodium and rubidium. Since these pioneer works all the alkali atoms 



V(R) « D + — + — + .. . 



(2) 




(3) 
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(Li, Na, K, Rb and Cs) (for a review see HiaiM!) then hydrogen [25], 
metastable helium |26j . calcium j2jj and ytterbium |2H] have been 
photoassociated. Preliminary results for heteronuclear alkali systems 
have also been reported [29} 130]. 

In a dilute medium, as the one present in the magneto-optical 
trap, the probability to find two atoms at a distance R is propor- 
tional to 4irR 2 e~ v ( R ^ kBT . Consequently as PA is a collisional process 
it is efficient only at large interatomic distance R. Therefore PA is 
particularly well adapted to the study of long-range molecular states. 

Because of the extremely narrow continuum energy distribution 
(on the order of ksT), the photoassociation free-bound transition be- 
tween the two free cold atoms (T ~ 100 /xK) and the ro- vibrational 
excited states is resolved at the MHz range (fcgT ~ h x 2 MHz). This 
leads to an extremely precise spectroscopy |24| . The kHz range has 
been achieved in rubidium starting with an atomic Bose-Einstein con- 
densate |31j . New available precise data from PA spectroscopy have 
stimulated the theoretical determination of more precisely values for 
the asymptotic coefficients jSH ESI EH • 

We shall present here new useful simple analytical formulas to 
extract the leading coefficient C n of the multipolar expansion within 
a 1% precision. To illustrate the importance of such a calculation, let 
us mention that this term occurs in the expression of atomic lifetime 
t of the first excited p level of a dialkaly molecule : 



3?ic 3 



4|C 3 k 



t (4) 

at 



where hui a t is the energy difference between the p excited atomic state 
and the s ground state. Indeed, a precise r value was obtained using 
a pure long-range expansion of the 0~(s + P3/2) potential curve of 
dialkalis [351 H3 133 EE] • 

This article is organized as follows. Section [21 is devoted to the 
fully detailed derivation of a first improved LeRoy-Bernstein formula 
using three new estimations respectively for the asymptotic part, for 
the repulsive branch part and for the intermediate part of the poten- 
tial curve V{R). In section El we take into account the next multipole 
coefficient C m . We finally obtain our main results, the general for- 
mula (|20|) for all the semiclassical NDE expressions and the improved 
LeRoy-Bernstein formulas (|21|) and (|26|). We apply these results in 
section0]on the 0~(6s + 6p 3 / 2 ) state of the cesium dimer (where n = 3 
and m = 6). We will discuss in great detail in the appendix El how to 
derive formula (J2J for all cases, so that this theory can easily be ex- 
tended to other long-range states. Indeed one goal of this article is to 
give a self sufficient theoretical background helping people interested 
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in using our new simple analytical NDE formula in the interpretation 
of photoassociation data. 



2 Improved LeRoy-Bernstein theory 

One of the simplest way to assign a given spectrum with a molecular 
potential curve is to isolate its vibrational progression and to extract 
an experimental C n coefficient, and then compare it to the theoreti- 
cal C n coefficient. This popular method makes use of the analytical 
semi-classical LeRoy-Bernstein formula (jSJ), that we propose here to 
improve. 



2.1 BKW assumption 

We use the Jeffreys, Brillouin, Kramers and Wentzel ((J.)B.K.W) 
semi-classical method and the Bohr quantization condition (e.g. see |39j ) 
for the vibrational level v at energy E = E v of a reduced mass \x par- 
ticle moving in a potential V(R): 

1 v^7 f" 



v + 2 = iJrJ R y (E - v(R))iR < 5 > 

R-(E V ) and R + (E V ) are respectively the inner and outer classical 
turning point of the vibrational motion (V(R-) = V(R+) = E). At 
the dissociation limit E = D the non-integer vibrational number v 
results of the formula is noted vd- 

For levels very close to the dissociation limit, the quantization 
condition is still a controversial subject |40 [ 1411 W2\ 143 j. For instance, 
it has been shown [11] that Bohr quantification condition should be 
modify at the dissociation limit by adding a 2 {n-2) term at the v + 1/2 
one. But the modification is of noticeable importance only for the few 
last levels (typically within less than 10 GHz energy range from the 
dissociation limit) of the potential where relativistic retardation 
effects or hyperfine structure appear, and where it is no more realistic 
to use the LeRoy Bernstein formula. Nevertheless, if needed we can 
furthermore improve the formula by taking into account this 2 (n-2) 
term or by using the third order semi-classical theory ^01 and adding 

d 2 V(R) 

a 487iV2/I (E-V(R))^ term in the quantization formula. 

However, as reported in the relative BKW accuracy is on the 
order of 

I t-iBKW Tjireali / irireal , 
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Thus, for levels close to the dissociation limit where typically v = 100, 
there is no need to improve the usual Bohr quantization condition 
to reach the 1% accuracy we are looking for. Consequently, in the 
following, we shall use the usual Bohr quantification condition @ 
and we shall see that other assumptions are less accurate than this 
one. 



2.2 Role of the non asymptotic part 

We define a "cut-off" outer-turning point where the potential V 
can be written as: 

R>R1 c 

V (R) « D + Jf ( 6 ) 

within a given precision. The potential V(R) and its asymptotic limit 
are represented in figure ^ Our goal is to reach a 1% precision then, 
if needed, R c + could be defined as: 

0.01 J^i- = -J^L (7) 
{R%) n {R%) m K ' 

With typical values as n = 3,m = 6, | C 3 1 = We 2 ciQ and \Cq\ = 
10000 e 2 ao, we obtain R c + 45 ao (where ao ~ 5.29 x 10 -11 m and 
e 2 = q l/{^e )). 

It is now possible to separate the non-asymptotic part from the 
asymptotic part (R > R+)- Taking the derivative of expression (J5J), 
we obtain (we use x = R/R + ): 

dv J2ii (-Cn) 1 ' 11 f 1 x~ n l 2 , 

= dx + 



dE 2-Kh (£) _ JR C + /R+ v 7 ! 



x 



1 **=±. (8) 



2irh Jr_ y/E - V{R) hto 

" v ' 

I dii \ n - a - 
KdEJ 

where the subscript n.a. is for non asymptotic and u is the classical 
pulsation of the vibrational motion. 

A physical insight on the role of the non asymptotic part can be 
obtained considering the classical definitions of the velocity v and the 
impulsion p: 

v= p _ y/2»(E-V(R)) 

Equation (jSJ can then be written as 4^ oc K We can see in figure 
^that the motion time is largely dominated by the asymptotic part of 
the potential, given by a multipole expansion as in formula (jHJ). This 
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distance A* 



Figure 1: The potential V(R) and the asymptotic potential D + C n /R n . Both 
potentials coincide when R is larger than the "cut-off" R c + . The repulsive 
branch is model by a linear formula for R < R°_. 

/ , \ n.a. , 

classical discussion tells us that ( J <C for levels close to the 
dissociation limit. The next step is then to restrict ourselves to the 
levels close to the dissociation limit 

2.3 Near the dissociation limit 

2.3.1 Correction in the asymptotic part 

For levels close to the dissociation limit we have the following inegality: 
y = R c + /R + -C 1. We can then write: 



Jy 



i x n / 2 x n / 2 f y x n / 2 

; dx = , dx — , dx (9) 

VI - x n Jo \fl -x n Jo \fl -x n 

= 0F r(l/2 + l/n) 2 ^ n + 2 

n r(l + l/n) n + 2 y v 6n + 4 y ; 

^ n + 2 

^ V? r(l/2 + l/n) 2 ^/ -C TO \- — 

n r(l + l/n) n + 2 1 +> \D — E ) 

where the main term (the only one taken into account in the "usual" 
LeRoy Bernstein law derivation) and the first correction term has 
been kept. This correction term brings a real improvement. Indeed, 
if in the last formula only the main term is kept, we need to take 
R+ = hR c + (y = 1/5) to obtain the integral Q value at a 1% level. 
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Concequently, with ~ 45 do, reaching a 1% accuracy level with 
the the "usual" LeRoy-Bernstein formula requires to use levels with 
R+ > 200 ao where the Bohr quantization problems occur. On the 
contrary, when using both terms, taking R + = 1.6i?+ is enough to 
reach the same precision level. 



2.3.2 Repulsive branch 

To express the non-asymptotic part in formula (jHJ), we will model the 
inner wall by a linear function using another cut-off R°_ (as indicated 
in figure and two parameters D and C: 

R<R C _ „ 

V(R) f» D + C{R C _ - R) (10) 

Other models (e.g. a potential with a 1/R behavior) do also lead to 
analytical formulas. 

The non-asymptotic integral in formula (jHJ can then be splitted in 

R c R c R c 

two integrals, using J R ~ = J R ~ + J R c + . The first integral is computed 
analytically using formula 1|1U|): 



1 :d^ = 2 ^^^^g (11 ) 



r_ y/E - V{R) C C 

We have moreover use the approximation E — D ~ D — D, because 
we are dealing with levels close to the dissociation limit (see figure 
Better accuracy could be achieved by keeping E in expression (jllj) . 



2.3.3 Intermediate region 

For the second integral J RC + , we use the assumption E — V(R) ~ 

D — V{R) that is valid in the intermediate region R c _ < R < R ^ 
(see figure Thus, this second integral becomes simply a number 
and does not vary with E. If needed, it can be computed using for 
instance a model potential like Morse, Lennard-Jones or a quadratic 
one. Finally the non asymptotic part becomes: 



dv 
~dE 



2irh 



2V D - D 
C 



+ 



Hi 



Rt VD - V(R) 



dR 



(12) 
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2.4 Improved LeRoy-Bernstein formula 



Using formula Q in equation (jBJ) and integrating it using expres- 
sion (|12j) leads to the improved formula: 



v D - v ps H~ 1 {D - E) V + 7 (D - £) 



(13) 



where 



(14) 



tt ft(n-2) r(2±l) 
and 7 is an extra parameter strongly related to R c + defined as follows: 
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2// 



(ra + 2)7r/i 
277 



2vr?i 



C 



+ 



(15) 



Ki y/D - V(R) 



dR 



Formula (|13() is a very simple one because all the three terms in 
formula (j!5j) match the same D — E behavior in formula (j!3|) . 

As vd — v does not depend on i?^, expression (|T1?)) shows that 
7 should also be independent of R c + . As a consequence (see expres- 
sion (EIJ)), the non asymptotic part (from R_ to i?^_) should follow a 

Tl_|_2 

behavior. This is not fully satisfactory as the non asymptotic 
part should not depend on the value of n which is a purely asymptotic 
parameter. This kind of trouble occurs whenever a cut-off is present 
in any theory. To prevent this caveat, we have hidden the cut-off R+ 
inside the only parameter 7. Thus, the final formula ()13|) is no more 
depending on the cut-off value R c + . 

Considering the added term ^{D — E) as a perturbation, it is pos- 
sible to reverse the formula l|13|) leading to our first improved formula: 



D-E 



VD 



2n 
V \ n-3 



Hn 1 



2n 



n 



VD 



V \ n-2 



H n 1 



(16) 



The LeRoy-Bernstein formula is then improved by simply adding one 
term depending on a single coefficient 7 that can be used as a parame- 
ter in a fit procedure. Will see in section |1] how this formula improves 
the fit to the experimental energies values. 



S 



3 Other multipole expansion coefficients 
General formulas. 



3.1 General NDE formulas 

To improve further the accuracy of the LeRoy-Bernstein formula, we 
can add other multipole expansion coefficients as in formula (j2j). The 
cut-off is then redefined so as to obtain with a typical 10% accu- 
racy: 

R>R1 C C 
V(R) ^ + D + ^ + ^ (17) 

To be more general, let us notice with LeRoy |l()j that the vibra- 
tional progression in vd — v, the rotational constant B v , the kinetic 
energy T, other BKW expressions of higher order, or higher rotational 
constants as D v , can be derived from the integrals ij^j: 



dR i J 



l(E) ~ L {E) #(E-V(R))W* dR (18) 



r R+(E) ( — 

JR-(E) R l {E - 

The derivative of relative to E is equal to —{k + 1/2)I % ^ 3 +1 v Thus, 

the NDE expression can be calculated only for the k = case: 
The method goes as follows. 

As a first step, we separate the integral in three parts J R ~ + f R <t + 
The first integral is analytically calculated using the linear expres- 
sion JTUJ) for V. Using l/y/E-V{R) « 1/y/D - V{R) in the inter- 
mediate region, we compute the second integral as a simple number, 
independent of E. If we use the same assumption as in formula 1)11)1 
the first two integrals, symbolically written as If' a '(E), can be approx- 
imate by the number If'^^D). In the asymptotic region (third inte- 
gral), where V is given by the polynomial multipole development (|17|) . 
the numerator in the wanted expression (|18j) is just a R polynomial 
expression, so it simplifies with the R l in the denominator. Finally, 
we need to calculate only a single expression: 

f R+(E) 1 

If{E)= = dR (19) 

lK ' Jb? + R l ^E-V{R) V ' 

the subscript " a " is for "asymptotic" and we use similar notations 
as |15j . We will only consider the case < I < 2. 

We simply have to follow the same kind of modification used in 
the previous section to calculate If{E). The computation, detailed 

in appendix ^ is based on first order correction in a c = 
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Consequently, using the notations j3 = ra+ 2~ 2 ^ anc ^ ^ = P~ i which 
can be null for the set of values (n, m, I) = (4, 5, 2) or (6, 10, 0) for 
instance, the complete (non-asymptotic) formula (|40|) can be written 

in a more useful and compact form L(E) = fn + !'^) — — , 1 dR: 

r IV > JR^(E) R i ^Je-V(R) 



( r< ^-1/2 r> /n p\ -& I P B {Pi 2) if 5 < 

<-G n ) C m j i HD _ E) - d5 = 

{5-±)B(6,D if<5>0 



where 75 groups all the constants terms as the non-asymptotic ones 
and depends on C n ,C m , R° + ; for a precise value, usually not needed, 
see equation (|40|). B(a,b) = p£~p^ are Euler Beta functions. 

Ignoring the 7$ term and choosing C m = 0, we recover the usual 
("non" -improved) NDE formulas. 



3.2 Vibration 

We shall first develop an example of vibrational progression v = 
— V^; ^ is nothing else than formula ©. The NDE de- 
velopment for vd — v is given by integrating the differential equation: 
iE = S'o, leading to: 



VD 



2 ^ r ( _ Cn) /3-l/ 2 £(Al/2) (D-E) 1 "? 



2tt% 



n 1 -0 



+j 5 (D - E) (21) 

(3B((3, §) if 5 < 

i ln(D - £) if «J = 
if 5 > 



/ ^< \<5-i/2 / n p\i-<s I P B (P'\) if 5 < 

+ {-Cn) ^ D ~ E l \\n{D-E) if 6 = 

n C n l — o 1 - 



where, for instance, 75 = -57^75 for 5 7^ 0. 



3.3 Rotation 

Using the expression of the averaged semi-classical probability of pres- 
ence (e.g. see 139]) for the radial wavefunction ib v : — , ^ M , and 

expression for uj using formula Q, we derive the NDE analytical ex- 
pression for the rotational constant B v = j^(if>v\jp\'4'v) [13] : 

B - ^ If < 22 > 

10 



It is well known that this formula is less accurate than the formula © 
for the vibrational progression. Indeed, we had to neglect in the non- 
asymptotic region (i.e. for small R values) the 1/R 2 term in the B v 
calculation. Our improved formula 1)20 [I should also help to solve this 
point. 



3.4 Kinetic energy 

Similarly, we can compute the average kinetic energy (T) = (ip v \E — 

(T) = lfr d +1 / 2 > = J^JLlll = (23) 

This problem is indeed crucial in photoassociation. During its de- 
excitation, the kinetic energy of the PA excited molecule is transferred 
to the two free atoms. They can then leave the magneto-optical trap, 
if their speed is sufficient, leading to a detectable signal. 



4 Experimental comparison 

The validity of the usual LeRoy-Bernstein formula (|16|) with 7 = 
has been studied ^3 but not extensively. We are going to study 
the own performances of our formula. 



4.1 Multipolar development for dialkalis molecules 

We shall not present here a complete review of the former applications 
of the LeRoy-Berstein law to photoassociation spectroscopy. Although 
we won't give an overview of the theoretical study of the Hund case 
(c) long-range potential curve, we want to give a brief introduction to 
the subject in order to be able to compare the calculated coefficients 
C n and C m with the modified NDE formulas. A detail introduction to 
the Hund case (c) potential curve calculation is given in appendix IBl 
Our final goal is to obtain the C n and C m leadings coefficients for the 
long-range states and to be able to take into account all terms needed 
to reach 1% accuracy in the C n value. We will focus on this section on 
alkali homonuclear molecules dissociating toward ns + n'p asymptote. 

Hund case (a) potential curve dissociating toward the ns + n'p 
asymptotic limit of the two identical alkali atoms leads to the following 
multipolar expansion: 



(-1 S~12^s 
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c 3 


9.997(23) 


c$ 


-17390 


/nil 

°6 


-11830 


U 8 


-5040000 


ri^s 


-16560000 


U 8 


-2256000 


Wis 

U 8 


-913100 



Table 1: Multipolar coefficient in formula (|24|1 for the cesium case at the 
6s + 6p asymptote. C3 is closely related to the dipole matrix element (J57|) 
and therefore the atomic lifetime (formula (@J)). The C3 value given here is 
extracted from atomic lifetime measurement Cg and Cg coefficients were 
theoretically calculated 



T/ — T/ — 3 _L 6 , °8 , 

where coefficients are given in the tabled for the cesium case. 

The Hund case (c) potential curve (see figure EJ) are obtained after 
diagonalization of the matrix M + 5M. M is given in table El SM is 
a correction matrix given in table El which has to my knowledge never 
been published before. After the diagonalization of such matrices we 
calculate the multipolar expansion, i.e. the power series by respect 
to 1/R. Results are summarized in table [2 Let us mention that the 
"real" accurate expansion will contain other terms coming from re- 
tardation (1/R dependence), Coriolis (1/R 2 dependence) or spin-spin 
(1/R 3 dependence) effects. But these effects are negligible compared 
to the C n /R n term in the multipolar expansion for the range of inter- 
nuclear distance (20 — 200 ao) we are working with. Nevertheless their 
contributions are evaluated in appendix [B] with the new relativistic 
lifetime correction 5M. 



4.2 Testing the improved LeRoy-Bernstein for- 
mula with the CS2 0~(6s + 6^ 3 / 2 ) state. 

Our group( |1B1 |H2 0231 ) obtained experimental photoassociation spec- 
tra with an accuracy of 150 MHz for all the seven allowed states (see 
figure El • We will focus here on the spectrum, from v = (D — E v ~ 
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V±2 U 
V±2 



V- 



±lu 



v± la = 



A 
A 



C 3 



^n 3 



_| _| 2_ 

i? 3 # 6 



2 



2+V7) Cj 9 A (2 (-1+y^) Cf+(2+V7) Cf)+2 (-10+7 yg) C 3 2 (1+e) 2 



-A 



3i? 3 
2C 3 (l+e) 2 
3R 3 



+ 



27V7AR 6 
9A(2C 6 n +Cf )-28C 3 2 (l+e) 2 

27. \ /."• 



A _ (-2+y^) C 3 9 A (2 (c"-Cf )+y^(2C n +Cff) )+2 (10+7n/7) C 3 2 (1+Q 2 
2 3R 3 + 27 v / 7Aij6 

A (-2W7) C 3 9A (2 (C 6 n -C|)+y^(2C 6 n +C 6 s ))+2 ( 10+7^7) C3 2 (1+e) 2 



„ + 

-A - 

A 
2 



3i? 3 
2C 3 (l+e) 2 



+ 



+ 



27 y/7 AR 6 
9A(2C^+Cf)-28C 3 2 (1+e) 2 



3 I ~" ' " 27 A i? 6 

2+v^) C 3 9 A (2 (-l+v 7 ?) Cff + ^+y^) Cf)+2 (-10+7^) C3 2 (1+e) 2 
3i? 3 + 27v / 7Ai? 6 



A _ 4C 3 (l+e) 
_/1 3R 3 



9A(2C 6 n +C 6 E )-4C 3 2 (l+e) 2 



A 
2 



5C 3 
3R 3 



+ 



27.1 /.'• 

9 A (cf+2 Cf )+4 C* 3 2 (1+e) 2 



, 4C 3 (l+e) 
_/1 "I 3R 3- 



27AR 6 



9A f2Cn+Cf )-4C 3 2 (1+e) 2 



27 A iJ 6 

5C 3 j_ 9^(g 6 n +2C 6 s )+4C 3 2 (l+e) 2 
27Ai? 6 



2 T 3R 3 T 



, A(2C 6 n +C b s )-4C 3 2 (l+e) 2 2C s n HC a S ' 
_/1 + 3Ai? 6 + 3R 8 



A I C 3 I 4 C "+ 2C 6 S )+4C3 2 (l+e) ; 
2 ~r fi 3 "t" 3Ai? 6 " t " 

A(2C"+Cf )-4C 3 2 (l+e) 2 2Cf^ 

^ "T Q /I E>6 I Q 



C,P s +2Cf s 



3R 8 



A 
2 



R 3 



+ 



3AR6 

A(c n +2Cff )+4C 3 2 (l+e) 2 
3Ai? 6 



3R 8 



+ 



C^ a +2Cf a 
3R 8 



Table 2: Two main term in the multipolar expansion in 1/R for Hund case 
(c) potential curve. The dissociation is taken at ns + n'p limit. Then, levels 
starting with A/2 dissociate toward ns+n'p 3 / 2 , while the —A levels dissociate 
toward ns + n'pi/ 2 - 
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15 20 25 30 35 40 

R <a„> 



Figure 2: Hund case (c) potential curves for the 6s + 6pj dissociation limits 
in the cesium case. The 7 attractive allowed curves (broad solid lines) for a 
dipolar transition, as the photoassociation transition from Qs + 6s atoms, are 
shown with repulsive (tiny solid lines) and forbidden ones (dotted lines). The 
curves are issue of a matching between short-range ab-initio calculation [UJ 
and long range ones given in the formula (|2*3jl and the table H 

77 cm -1 ) to v = 133 (D — E v ~ 0.4 cm -1 ), of the external well of the 
0~(p3/ 2 ) state. To test our formulas we will use the data extracted 
from a RKR study and published in |51) . Using tableland [2 we shall 
take n = 3, m = 6, 1 = and C n —10 e 2 ao, C m 65000 e 2 ag. This is 
one of the most ruthless case to test a LeRoy-Bernstein type formula 
because almost all the assumptions used in its derivation are wrong 
or could be discussed. The improved version shall then be needed. 

Firstly 5 < (and C m > 0) lead to problems that LeRoy in 
reference ^S] had. This is due to the fact that Tq '^(0) = oo in formula 
(|38|) . Our derivation avoid them using T f(y n ) which is finite. 

Secondly n is small and we use a "pure long range state" where 
R is always large. Therefore, x n = (R/R + ) n is not so small and 
keeping only the first order, as done in the "usual" LeRoy-Bernstein 
derivation, in the series (JHJ) might be not accurate enough. 

Thirdly, the term C rn as 65000 e 2 ag is larger than in usual Hund 
case (a) potential curve, therefore the C m /R m correction term is also 
large. Furthermore is of greater value and y = R+/R+ is not small 
(see for instance formula ©)• 

Fourthly, the inner wall is smooth and less steep than usual (see 
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figure [2J) . Consequently a small phase is accumulated on the inner 
wall (see formula 

In the fifth place, the potential curve is only 80 cm -1 deep. So, 
1/y/E — V(R) is quite large and the non-asymptotic part, defined in 
expression (|H]>. is therefore important. For the same reason, the as- 
sumption: l/\/E — V(R) ~ 1/y/D — V{R) in the intermediate region 
can also be wrong. 

For R + rj 40 ao, i.e. D — E m 30 cm -1 , the next asymptotic 
coefficient C 8 /R s is already 4% of C 6 /R 6 (it reaches 10% for R « 
25 ao). Therefore a choice of i?+ 35 — 40 ao is probably good enough 
to obtain \a c \ < 0.1 and a non-asymptotic part not too large. As 
mentioned in section r2.3.11 we also need R+/R+ large enough to get 
a good precision in our fit. Concerning the percent accuracy goal, the 
discussion following formula Q has indicated a restriction for the fit 
at R + > 1.6i?^ 55 ao corresponding to D — E < 12 cm -1 . 

For a physical insight, we give in table El some typical values for 
all the terms present in formula (|21|) . The potential curved used to 
numerically evaluate all the terms in table 01 including the term 

a/2^ / 1/2 (R%) n(3 1 ^1/2 Cm (R+) US \ 

15 = ^h\ l1 / ^^ + 2 ( - c " ) c n ^r) 

(25) 

is simply the diagonalization of the g matrix describe above where 
the Cs and e coefficients are assumed to be zero. These results confirm 
the well known fact the "usual" LeRoy-Bernstein formula (with the 
sole term (D—E) 1 '^) won't gives results better than one percent. This 
also confirms that the "usual" LeRoy-Bernstein formula is by chance 
much more accurate than it should be because the non-asymptotical 
parts If' a '(D) in 75 almost perfectly canceled with the other terms 
such as 73. 

In the non-asymptotic part we have made the approximation, as 
in formula (Jr2|). If- a -{E) m If a -(D) which, for instance with D - E m 
10 cm -1 , leads to an accuracy of the vd — v value of about 0.3. This is 
of similar importance than other contributions listed in table As a 
consequence, it is useless take into account contributions smaller than 
0.3. Therefore, this means that in order to improve our formula, we 
shall not incorporate second order terms (such as 0(y n )) but we shall 
rather have to take into account more carefully the non-asymptotic 
part. This can not be done without adding other unknown parameters 
in the development as 4%l or i?^-dependent terms and without keeping 
E in formula Ijllj) . In a sense our theory is the best with only one single 
unknown parameter added 75 to V£),D,C n and C m . 

Our theory, see expression (|2*T|) yields vd — v = f(D — E), but we 
would prefer to adjust the theory to the experimental energies, i.e. to 



15 



D - ^(cm- 1 ) 


30 


10 


1 




172 


143 


97 


T- & iD) = &£If-*-(p) 

i \ / Z7T/1 L y / 


12.7 


4.2 


0.4 




16.2 


4.5 


0.4 


7/3 


-9.8 


-3.3 


-0.3 


is 


3.7 


1.2 


0.1 


(D - Ef- d 


-1.8 


-0.5 


-0.03 


0(y n ) 


-0.7 


-0.1 


-0.001 



Table 3: Estimation, for several D — E values (in cm -1 ), of the terms involved 
in vd — v formula (J2H)- We use the analytical Q~{pz/2) potential curve with 
C$ = e = 0, where we have chosen R c + = 35 ciq. Terms are labeled by the 
D — E power in the formula (|21|). The three separate terms in formula (|25|) 
are respectively noted 7 n - a -, 7^, 7$. We also evaluate the contribution from the 
neglected 0(y n ) term in equation (|36|) which is the second term in formula 



fit using D - E = - v). Table 01 indicates the formula {20 

is largely dominated by the first term H~ 1 (D — E) 1 ^^ and that we 
could safely, except for very large D — E values, use the only first 
order inversion procedure used for instance to derive formula (fTB|) to 
find an accurate enough general reversed formula: 



D-E 



VD 



1 

1-/3 



H n 1 

(_C7 n )*-i/2 C , 



[1 



1 



1 



1 



- (3v D 



!~ ( VD 
>75 



1 

1-/3 



+ 



n 



C n 1 — 5 V H n 1 



1-6 

1-/3 



V H n 1 

< iln(^?)^ 
k (5-1)5(5,1) 



This formula gives an explanation for the origin of the Pade coeffi- 
cients used in " usual" NDE formula [H2j > as long as an explanation for 
their values. Indeed, Pade formulas assume a mathematical, without 
a-priori a physical meaning, polynomial quotient expansion in vd — v 
for physical value as vibrational series or rotational constants |16j . 
Our formula leads directly to such polynomial expansion and gives 
physical interpretation for the coefficients. 

In table HJ we present fit results (done with Mathematica software, 
with 100 iterations in the non linear fitting procedure) for all the vibra- 
tional levels and for the vibrational —30, —10, —5 and— 2 cm" 1 of the 
0~ state. It should be remembered that we use data from RKR compu- 



(26) 



if 5 < 

if 5 = 
if 5 > 
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tation where C3(0„ ) was kept fixed at -10.47e 2 ag 0]. The accuracy 
is also much better than the experimental one because we use a NDE 
theory to fit the NDE-RKR calculated levels and not the experimental 
levels. The table 0] shows an improvement when using our formula, as 
opposed to the usual LeRoy-Bernstein one. Our method is much more 
stable than the usual ones when the fitted energy range changes. Thus 
we are able to extract C3 ~ —10.5, vd ~ 214.8, D ~ —0.09 cm" 1 which 
are very close to values found by a complete NDE-RKR analysis |48| 
C 3 = -10.47, v D « 214.6,1) « -0.09 cm" 1 . 

Our method seems suitable for extracting asymptotic coefficients 
at the percent accuracy level. The method could be accurate enough 
to need the corrections factors as the 5M matrix as the e one which 
has never been used up to know. The method can also be applied to 
other states as the 0~ (P1/2) or the l g {jp^/2) ones. The Ogipi/2) state 
has a 1/i? 6 asymptotic behavior (see table |2j) and then should gives 
information concerning the next multipolar coefficients Cq and Cq-. 
The l s (p3/2) leads to n = 3 and m = 8 in the cesium case because, 
see table 01 and value given in table Q the 1 /R 6 term is accidentally 
very small and therefore negligible comparing to the 1/R S one. 

5 Conclusion 

We have derived general improved NDE expansion formulas (|20|). in- 
cluding LeRoy-Bernstein one (|2fjjl. leading to a better accuracy in the 
determination of the asymptotic coefficients. 

Such expressions can be useful for further photoassociation exper- 
iments to extract the asymptotic coefficient C3 or Cq and the atomic 
lifetime. 

The method gives also a physical meaning for the Pade coeffi- 
cients used in usual NDE formula. The method could then be used 
as a starting guide for Pade coefficients calculation. Furthermore our 
theory includes, as physical parameters, the two leading multipolar 
coefficients C n and C m . The single added parameter 7,5 contains in- 
formation on the repulsive branch, the intermediate internuclear dis- 
tance behavior and analytical calculation of the non-asymptotic part 
of the vibrational phase in a C n /R n potential curve. We have shown 
it is not reasonable, without adding another parameters, to develop 
further than we did the approximation in series development of the 
analytically known asymptotical part. 

The author thank C. Amiot, H. Blanchard, O. Dulieu, D. Hardin 
and N. Vanhaecke for many helpful discussions. 

^Laboratoire Aime Cotton is associated with Universite Paris-Sud. 
web site: www.lac.u-psud.fr 
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method 


range 


D 


C 3 (07) 


vd 


75 


c 6 


a ht 


© 


all 


0.85 


-17.94 


245.1 








9879 




-30 < E v 


0.16 


-14.10 


228.7 








1906 




-10 


-0.03 


-12.16 


220.6 








230 




-5 


-0.06 


-11.60 


218.4 








59 




-2 


-0.08 


-11.12 


216.6 








11 


(HED 


all 


0.06 


-12.42 


222.6 


14686 





1272 




-30 


-0.06 


-11.27 


221.7 


20142 





151 




-10 


-0.084 


-10.79 


215.7 


26378 





14 




-5 


-0.088 


-10.68 


215.3 


29544 





9 




-2 


-0.091 


-10.51 


214.8 


37070 





9 


(BSD 


all 


-0.080 


-10.71 


215.6 


56034 


-1.310 7 


61 




-30 


-0.089 


-10.55 


215.0 


62565 


-1.510 7 


10 




-10 


-0.090 


-10.49 


214.8 


69581 


-1.710 7 


9 




-5 


-0.090 


-10.48 


214.7 


73770 


-1.910 7 


9 




-2 


-0.091 


-10.52 


214.8 


37019 


-25061 


10 



Table 4: Value of the coefficients D, C 3 (0~),Vd, Is, C 6 for the 0~(6s + 6p 3 / 2 ) 
cesium state, the numerical fit was done using formula (|26|) with k = 5 
parameters for different set of N vibrational states determined by the range 
of available E v values. When yg = Cq = (k = 3), we recover the usual 
LeRoy-Berstein formula (J3J). When C§ = (k = 4), we recover the first 

improved formula (HEJ). a fit = j^^Ef=i{E v - Ef) 2 is given in MHz units 
so as to be compared to the experimental accuracy o = 150 MHz. 
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A Derivation of the improved LeRoy- 
Bernstein formula 

We will detail derivation starting from formula (|19|) . 
Using similar notations as ^3], we define: 

x = R/R + 

y = R~ + 

a(R) = -^j^ = a c (y/xr- n where (27) 

a c = a{R c + ) and (28) 
a = a{R + ) = a c y m - n (29) 

C m /R m <C C n /R n in the asymptotic region (R > R+) is equivalent as 
saying that a c <l. This inequality can be assured by a right choice 
of R c + . We shall then use a development in Taylor series about a c = 0. 
We define the "zero" -order parameters: 

y °" = (-c n )/{R%r m 

in — n 

„ „,m-n C m (D-E\ n 
ao = a c y = — ^ _ I (31) 

Then expression (|17|) becomes: 

R>R1 C 
V (R) * D + j±(l + a{R)) 

and equation E = V{R + ) leads to the following implicit equation for 

y- 

y n (l + a c y m - n ) = yZ (32) 
Using (3 = n+ 2n 21 > ^> ec t ua ti° n (|19() is easily written as: 

/ p?c \n/3 r i n/2-l i 



The derivation is quite similar to that of equation As < x < 1. 



we have 1 < ^~ g n < — ■ As a consequence, the second term is a 
perturbative one (because a(R) < a c <C 1), so we can expand the 
inverse of the square root in its converging Taylor series: 

It = (-C n )-V 2 (^Yjti ~\ /2 ) (a c y m ~ n ) k T^ m {y n ) (34) 
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where ( ) = " 1 ^ 2 ^ ^'"^ k+1 \ using u = x n and: 



W = »t^r d " (35) 



We use a similar notation as in LeRoy's paper [T2]. However, contrary 
to his computation, where y = is fixed before expanding into Taylor 
series in a (not a c as we did), thus leading to diverging integrals 
T^ m (0), we isolate here the dependence in a constant term a c and 
take into account the fact that a is y dependent. Indeed formula (|29|) 
indicates a — > when y — > 0. Furthermore, LeRoy can not deal with 
a < (i.e. C m < 0), as we do in section [Ol 
Equation @ has shown that 



= ^ (A1/2) - ^(1 + 0(,/< ')). (3fi) 

n nfj 

with the Euler Beta functions: 

B(a,b) = I u a - l (l-u) h - l du ; a,b > 

JO 

r(o)r(6) a + 6 



r(a + 6) 



-5(0 + 1,6) (37) 



Similarly we have to expand the integral (JHSJ) into a series about y n = 
and then to calculate the integral from to 1. Inserting 1 = (1 — u)+u 
and integrating by parts, we obtain the first term of the asymptotic 
series about y n = 0: 

ryv)= 1 } { -iny" if 5 = 

{ (1 - 26)B(6, |) + 2pB(fi, |) if 5 > 

(38) 

with 

^ — m 3n — 2m — 2/ + 2 

n 2n 
Let us focus hereafter on the first order correction in a c . Equa- 
tion (1321) leads to y n yg (1 — ao) and the formula becomes: 



' 4M r^ m ( 2/o «(l-ao))(l + /3ao)-^o^ m ( yo "; 

(39) 

where we have to consider only terms in formulas (|36|) and 1)38(1 with- 
out the 0(y n ) term. We have neglected terms that would be part of 
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the correction in second order in a c , and would lead to an explicit 
dependence of the NDE formulas on the cut-off (which then won't 
be no longer hidden in the 75 term). But, this might not be the best 
strategy. We know that y = R+/R+ is not necessary small (but y n 
will be). It might be better to approximate the integrals (|33|) by Pade 
series in y n than by Taylor series. Another method would be to obtain 
an analytical solution for the integral using for instance a third 



-C 1 



-) in the multipole expansion Q17|). 



term (as for instance c Riri 

Finally, with those choices, the full (non asymptotic) formula reads: 



II 



.^-1/2^09,1/2) (D-E\- p 



n 



irP)-(-Cn)- 1/2 ^^) + 



(40) 



(_ gn )-i/2 Cm / D-E 



D-E 



if 5 < 

2 iiLy-L^ J-i ) 1 2 ^ —C if ^ = 

(i-i)fl(U)-W{) if 5 > 



1 I 28 \-C n /(Rir 



where to summarize our calculation of formula (fTH|) , we have used the 
assumption E ~ D in the non-asymptotic part J R + . 



B Long-range states. 

We shall give a brief introduction to the Hund's case (c) states and 
potential curves calculation for long-range molecules. We shall mainly 
discuss the case of neutral alkali-atoms. The discussion can easily be 
extended to all atoms. We refer to the following books and articles [331 

EH EH ESI EE1 EZl ESI ES] 

B.l Non relativistic electrostatic interaction 
and multipolar development 

In the following, we shall discuss the interaction between two atoms A 
and B (radius ta and r&) formed by a core A + and a valence electron 
ei for the first one and by a core B + and a valence electron ei for the 
second one. 

The first assumption for the Born-Oppenheimer Hund's case (a) 
potential curve calculation is to neglect all the relativistic parts. Then 
the atomic Hamiltonian for atom A leads to 1 1 : nlmim s ) a states 
(quantize on the internuclear axis) represented by the eigenfunction: 

^l m (r A i)Xs=l/2,m s = ^^Yr(f A i)Xs=i/2, ms (41) 
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where ?ai are the polar angle of tai- Then the eigensystem of the full 
electrostatic Hamiltonian 



H = Hf + H% + H cl (42) 

can be computed. We will use the first and second order perturbation 
theory for H el . 

The large internuclear distance assumption (R » ta + Tb) leads 
to a Taylor series about 1/R. It is convenient to define the 2 J polar 
irreducible tensors Q 3 by (for atom ^4): 



Q7 = ^^^AiYr(?Ai) (43) 

For a numerical computation of the matrix element, this expression 
should be multiply by 1 + e(rAi) in order to take into account the 
effect of all electrons and not only the valence one |32| . 

With these notations, the Taylor series is H el = Eij=o 

where 



in g) i - 1 )J(i + j )lQ T (A)Q- m (B) 



^{i + m)\{i - m)\(j + m)\(j - m)\ 

It is then straightforward to see, for instance, dipole-dipole interaction 
has a -R -3 behavior. 



B.2 Molecular symmetries and Hund's case 
(a) states 

It is useful to use the molecular symmetries for H . Because for us the 
nuclear spin is passive (see HHIED] for details), we will not study the 
symmetries for the total (electronic, rotation, vibration, nuclear spin) 
wavefunctions, as the permutation for bosons or fermions nucleus, but 
only for the electronical part: 

• Orbital electronical rotation around the molecular axis Oz, lead- 
ing to a well defined projection firing = HA on Oz. 

• Orbital electronical reflexion a' v versus yOz (or a different plan ^56 ) 
(eigenvalue a' = ±1). This commutes with the previous one if 
A = 0. 
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• In the case of identical nuclear charge = Zb- electronical or- 
bital inversion I el f(rAi,rB2) = /(— reii — ^42) with eigenvalues 
u' = ±1 and states noted <? (gerade) for lo' = 1 and u (ungerade) 
for a/ = — 1. 

• Rotation of the electronical spin S. So 5 2 and S z are eigen 
operators with h 2 S(S + 1) and Tims = ^E eigenvalues. 

The electronical spin is passive so it is easier to separate it and 
to define molecular spin states \S,ms = E) by: 

|0,0) 
|1,-1> 

|i, o> 
IM> 

with the convention m s is noted in first place in the \m s m' s ) 
notation. 

Spherical harmonics formulas |61| for the atomic wavefunction (|41|) 
lead to: 

cr v \l : nlmim s ) a|2 : n I mim s ) b = (45) 
|1 : nl — mim s )A <8> |2 : n Z — m\m^)B 

I el |l : nlmim s ) a\2 : n'Z'm^m^)^ = (46) 
(— l) r |2 : nYmJm^xIl : nlmim s ) B (47) 

Electrons are fermions so the final state has to be anti-symmetrical 
for the electronical exchange P12 = P\% P^ m : 

P12II : nlmim s )A ® |2 : n'l'm'im' s ) b = 
|2 : nlmim s )A (8> |1 : n' I' m'im' s ) b 

Finally, with the usual Hund's case (a) notations, the molecular 
state we want to calculate is: 

2S+l A (^') \ 

n,(u/) / 

where = A + E. Some other basis, as the Wang basis |56j with |A|, 
can also be used, as all calculations can be done in any complete basis 
; the one we choose leads to simple expressions. 

We are focusing on this paper to the states reached by photoassoci- 
ation. Most of the photoassociation experiments start with two atoms 
in ground state (ns + ns state) photoassociated toward the first ex- 
cited asymptotes ns + np. Therefore, we shall continue the discussion 
with one atom (A or B) in nl state and the other one in n'V state with 



= ^(l + ->-|-+» 
= ^(l + -> + |-+» 
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I or /' null (but n 7^ n' is not required). It is nevertheless simple to 
obtain the formulas for the general nl +n'V configuration |55 [ I60 [ 162] . 

From the previous expression for symmetry operators, it is easy to 
verify that the wanted expression is: 

2 ^A« (ns + n'l'))°= (48) 

1 -i- ( 1 \& P orb 

^ { ' 12 -c(\l:n00) A \2:n'l'A) B + (49) 



V2 

uj'(-1)- s ~ 1 '\1 : nl'h) A \2 : n00) B ) 



\s,n- A) 



where the exponent means we did not use the perturbation theory 
yet, so this state is the zero order state for a given internuclear distance 
R. c is a normalization constant slightly i?-dependent due to the 
exponential overlap between ^ A (f + R) and H> B (r). Similarly , the 
P°2 b operator will put the electron e\ (resp. e-i) close to the core B + 
(resp. A + ): this leads to another exponential correction, known as 
the exchange correction H21 EU • 

We will work with large enough internuclear distances to avoid the 

overlap and exchange terms. We can then make use of — -j= 12 = 1 
in expression (|49[) for the multipole coefficient calculation and c = 
l/\/2 when (n,l) 7^ (n',l f ), c = 1/2 otherwise. Several authors tried 
to estimated what is the right cut-off to safely neglect the overlap and 
exchange effects. Let us just mentioned the article |65j modifying the 
simple 1973 LeRoy limit: 

R > 2(r A + TB ) 

leading to typical values close to 20 do • 

Among the 24 ns + n'p states only 16 will be non-degenerate as 
shown in figure 121 

B.3 Interactions, exchange and Hund's case 
(c) curve 

B.3.1 Interactions and multipole coefficients for Hund's 
case (a) curves 

The next step, starting with this zero order basis, is to apply the 
perturbation theory, or better, the degenerate perturbation theory 
(see [HH1II2]), to the H el perturbation. 

Calculation is straightforward using formula l|44|). As before we 
detail only formulas for the ns + n'p asymptote. This choice yields 
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for the first order perturbation to i = j = 1, i.e. a dipole-dipole 
interaction, and energy: 



V„ 



2S+1 a(+) 



(50) 



/2S+1 a(+) 



rel 



if 

j 1 + |A| i? 3 



2S+1 a (+) 



A nL' ns + n P 



where 



C 3 



, and the dipole is 



\l2<Ke ) 

(ns || Qi || n'p) = — (n'p \\ Q\ \\ ns) 
V3(n'pmi = 0\q e z\nsmi = 0) 



(51) 



Second order perturbation theory leads to the so called polarization 
terms (London or dispersion, and Debye or induction) yielding the 
final multipole expansion formula: 



V„ 



C' c f 



(52) 



where E p is the dissociation limit energy ns + n'p. Theoretical value 
for several C 6 and C' 8 coefficients can be found in . 



B.3.2 Hund's case (c) versus Hund's case (a) states 

For heavy atoms we must add another perturbation term, the spin- 
orbit one, in the Hamiltonian ^56 j 

^so = A so lt^+^4 + hzSu + 1 ^ 2 

where A so is the atomic spin-orbit constant. Due to spin-other orbit 
{e<i on A + and e\ on B + ) or mixing with curves coming from n"l" + 
n"'l"' dissociation limits, A so is in fact slightly R dependent. 

The new Hamiltonian results in less symmetries. Only J = L + S 
is a good operator leading to fl as a good quantum number. Fur- 
thermore the electronical reflexion a v has to act on spin also. It is 
straightforward jHSl IMj to see its eigenvalues verify a = (—l) s a'. It 
is then better to work with a new basis: 
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Table 5: Matrix given fig u Hund's case (c) eigensystem (R variation is not 
written) for ns + n'p limit. With A = h 2 A so = (E Pa/2 - E Pl/2 )2/3. 



± 2 



^( 3 n ±2 , g , u ) + 1 



±1 



Ogtu 



/ v( 3 n ±lig , u ) 

\ 2 2 



±A 

2 



V ('''>- : I. ' 



n 3n o,g^ ' 



-A 
V2 



2 V2 



^( 3 no, g ,u) 2 %. 



A 
V2 



This definition is valid for {A, Q} / {0, 0)}, on the contrary 2S+1 ^qI^ 

is already an eigenstate for a = (— l) s . 

Using H so definition and formula (|49|) yields the block matrix lxl 
for |0| = 2, 2 x 2 for |0| = and 3 x 3 for |fi| = 1. Finally, we only 
have to diagonalize these matrices, given in table to have Hund's 
case (c) states fijy and potential curve versus the known Hund's case 
(a) ones. 



B.3.3 Hund's case (c) versus Hund's case (e) states 

In heavy alkali atoms, H so is quite larger than H cl at large inter- 
nuclear distance. Consequently it is better to work with a "fine 
structure" Hund's case (e) basis |1 : nljmj}^ <8> |2 : n'l'jm'^B where 
Hf + Hz + H so is diagonal. 

To calculate the only missing perturbation H , we need, as before, 
to find a well (molecular)-symmetrized basis |34| I55j . To avoid such a 
work, we take advantage on the already known matrix given in tablc[SJ 
By definition the Hund's case (e) basis is just formed by the eigenstates 
of the matrices given in table[3when only the spin-orbit is present (all 
the electrostatic interactions Vn, Vs are set to zero). 

Let us just illustrate it on the 0~ case. The transition matrix is: 
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Then the desired matrix of H cl in the Hund's case (e) basis is just 
/if°9 _ PM° a p- 1 

element related to n'pj [6*T] : 



= PMy H so +H e1 ^' 1- ^ * s then better to use reduce matrix 



D3/2 = {nsy 2 || Qi || np 3 / 2 } = -^\ ns II <2l II ( 53 ) 

a/2 

= ~/^( ns II II n '^) ( 54 ) 




1 ^3/2 1 ^3/2^1/2 

M °s = [ ^3/2 ~ D 16^o 8tt £0 _ & ] (55) 



E 



where, for sake of simplicity, we have kept only the C^/R s leading 
term in expression (|52j) . As expected the spin-orbit is diagonal leading 
to the dissociation limit toward ns + n'p 3 / 2 with an energy E p3/2 = 

E p + h and toward ns+n'pi/2 with an energy E Pl/2 = E p — h 2 A so . 



B.3.4 Relativistic correction 

This new matrix is also useful, as we are going to see, to add a rela- 
tivistic correction to the C3 coefficient in the case n = n' . 

Indeed, the experimental lifetime measurements noted r 3 / 2 and 
T1/2 for npz/2 and npi/ 2 jlfi j Ifi7| lfi%] disagree with the predicted ratio 
-D3/2/-D1/2 = V2 (see formulas (|53]l and (jHljl ). This relativistic cor- 
rection has to be taken into account. Let us define the e correction 
by 



B 3/2 'V2 \| 2t 1/2 £ 1 2 /2 V 1.9800(0) 

where the experimental value is given for cesium. A new definition for 

C 3 = ^^\(ns 1/2 \\ Qi \\ np 3/2 )\ 2 (57) 

leads to the same value than the previous one (see formula (|51|l ) for 
e = (the value for cesium is C 3 = 9.997(23) gS]). 
Then the new matrix (|55|) is: 
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Using the transition matrix P we can then include the perturbation 
in the Hund's case (a) matrix to find: 

V — e # + ^r- ^ P -2 — j^h- / 

Comparing this matrix to the one calculated without the e correction 
leads to the correction matrix 8Mr a ) given in table H3 In the cesium 
case e ~ 0.005 therefore the corrections factor (1 + e) (or even more 
(1 + e) 2 ) are (just) needed to have the percent accuracy we are dealing 
with. 

Their is a second relativistic correction known as the retardation 
effect [§S EH EE] ( see also reviews (721 E3 and recent articles [TUTS] ). 
The main retardation correction for E (respectively II) states concern 
the C3 coefficient which has to be multiply by /s (respectively /n) 
where: 

R R , R 1 R 2 

= COS I + I Sm l" 1+ 2A^ 

R R R R R 1 R 

fu = -^^I + I^X + ^I- 1 -^ 

where A = c/to (2400 ao for Cs) and Tiuj = E p — E s . This theory is 
limited for R < ct (several centimeters) due to the photon lifetime. 

We will consider this correction as almost negligible for our purpose 
of 1% accuracy because our description will focus on R < 200 ao ~ 
A/12. 

To conclude this calculation, and in order to obtain a precise poten- 
tial curve determination, we have to include some other small effects. 
These effects are usually negligible to obtain the two leading terms 
in the multipolar extension as it is needed in our NDE expressions. 
Therefore they should not have any incidence in our calculation. 

1. Spin-orbite R (dependence), exchange, overlap. 

These effects will mainly affect the intermediate part of the po- 
tential curve and not the " pure-" long range part we are interested 
for our asymptotic calculation. 

2. Spin-spin. 

As previously discussed, the spin relativistic effect has to be 
taken into-account for a precise potential curve determination |71j . 
Spin-rotation, dipole-(spin-dipole) are negligible. The spin-spin 
interaction leads for instance for the 0~ matrix to the correction: 

h 2 e 2 , . w , 

(59) 
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Table 6: The relativistic correction SM to the matrix proposed in table (|5|); 
e value is defined by formula States g correspond to a;' = 1 and w to 

c/ = -i. 



5±2 



<5 ± lg,u 

50 



g-u 



50 



g.u 



\Ul'+l 



9R 3 
3.R 3 



/ 2(e-3) 

±2e 
V-(3 + 26) 

-4(3 + 2e) 
-^(9 + 46) 
-4 -v^ 
-V2 4 



±2e 
2(3 + e) 
T(9 + 2e) 
-V2(9 



-(3 + 2e)' 
T(9 + 2e) 
2(6 + e) 

-4e) 



-4(3 + e) 



This is a negligible term in the multipolar development because, 
for cesium, ss (7 3 x 5.10~ 6 . 

3. Rotation and Coriolis. The rotational part is given by 

rr rot _ _ (J-L-Sf 

2fiR 2 2fiR 2 

This usual derivation j^E] yields for instance in the 0~ case the 
matrix correction: 

h 2 < J(J + 1) + 2 2^2 \ 
2ixR 2 \ 2^2 4 + J(J + 1) ) 

For R ~ 200 ao we have 2 ^ Ri ~ 10~ 4 ^|. This correction is also 
negligible for the multipole expansion. 

4. Kinetics coupling and mass polarization terms. 

These terms lead typically to a correction of less than one per- 
cent [Tfij . 
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